Instability Regions in the Upper HR Diagram 
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Abstract 

The following instability regions for blueward evolving supergiants are 
outlined and compared: (1) Areas in the Hertzsprung- Russell (HR) dia- 
gram where stars are dynamically unstable. (2) Areas where the effective 
acceleration in the upper part erf the photospheres is negative, hence di- 
rected outward. (3) Areas where the sonic points of the stellar wind 
(where v win d = v 9aU nd) are situated inside the photospheres, at a level 
deeper than tr om * — 0.01. We compare the results with the positions 
of actual stars in the HR diagram and we find evidence that the recent 
strong contraction of the yellow hypergiant HR8752 was initiated in a 
period during which < g e /f > < 0, whereupon the star became dynami- 
cally unstable. The instability and extreme shells around IRC+ 10420 are 
suggested to be related to three factors: < g c ff > < 0; the sonic point is 
situated inside the photosphere; and the star is dynamically unstable. 

Key words: stars: atmospheres; interiors; supergiants; evolution; 
HR-diagram 
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1 Introduction; hypergiant instability 

The apparent instability of many stars in the upper part of the HR diagram 
has different causes, dependent on the stellar properties, which, in turn, are 
partly related to their locations in the HR diagram. Observations obtained in 
recent years are indicative of various modes of interior or atmospheric insta- 
bilities among yellow hypergiants and S Dor stars (=LBV’s). Evidences for 
these instabilities in yellow hypergiants were summarized by de Jager (1998). 
Specifically, they refer to phenomena such as the enormous pulsational ampli- 
tude (A R ~ 0.25 R) of p Cas (Lobel et a!., 1997), the ’bouncing’ evolutionary 
motions of p Cas (reviewed by de Jager, 1998), of IRC-l-10420 (Oudmaijer et 
al., 1994), of Var A in M33(Humphreys, 1975, 1978) and of HR8752 (de Jager 
and Nieuwenhuijzen, 1997). Another indication of atmospheric instability may 
be seen in the extended clouds of gas and dust around IRC-h 10420 (Mutel et 
al., 1979; Jones et al, 1993; Oudmaijer et al., 1996; Blocker et al., 1999). In 
S Dor stars the large outbursts and the quasi-oecillatory temperature changes 
(cf. review by van Genderen, 2001) have been suggested to be due to their dy- 
namic instability (Stothers and Chin, 1996). The understanding of the causes 
of these instabilities may profit from a delineation of areas where one or the 
other mode of instability prevails. It should be noted in this connection that, 
when speaking of stellar instability, we do not refer to the observed microvaria- 
tions (reviewed by van Genderen, 2001) whose origin might possibly be sought 
in the excitation of strange mode pulsations or other low-amplitude radial or 
non-radial oscillations. 

We summarise earlier work in this field. 

Stothers and Chin (1996) showed that in certain areas of the HR diagram 
the mean value of T\ (= (dlnP/dlnp) a rf) can take values below 4/3, which implies 
dynamic instability of the star. As a consequence, a highly evolved star to which 
this applies can be triggered to a phase of steady expansion or contraction. Thus, 
they were able to define two regions of dynamic instability in the upper part of 
the HR diagram, one for Teff< 10000, another for higher temperatures. These 
areas were called the ’yellow-red’ and the ’blue’ dynamic instability regions. 

Nieuwenhuijzen and de Jager(1995, summarized by de Jager, 1998) outlined 
two regions in the HR diagram where in the atmospheres of blueward evolving 
stars, hence very evolved objects, five conditions are obeyed. These are: g e fj < 
0.3 cms -2 ; dp/dz < 0 (z is the vertical ordinate) in the relatively deep parts 
of the photospheres; the sonic point, Le. the level where v win d — v a0 und > lies 
inside the photosphere; the sum g e jf 4- g^u < 0 during part of the pulsation; 
and Ti < 4/3 in part of the line-forming part of the photosphere. Their two 
regions were baptised the ’yellow void’ and the ’blue instability region 1 . 

In our studies of supergiant instabilities it became clear to us that it may 
advance the understanding of stellar instability (or quasi-instability) by consid- 
ering the various causes separately, lb that end we will delineate the regions 
in the HR diagram where supergiants or their atmospheres are unstable in one 
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way or the other. 


2 Areas of stellar dynamic instability 

Ritter (1879) showed that for radial dynamic stability the ratio of specific heats 7 
should exceed the value 4/3. More generally, Ledoux (1958) found that for a real 
star the first generalised adiabatic exponent T 1 , suitably averaged, should exceed 
4/3 in order that the star be dynamically stable. Following that line, Stothers 
and Chin (1995) and Stothers (1999) demonstrated that a non-adiabatic, spher- 
ically symmetric envelope of a star is dynamically unstable when o 2 < 0, where 
o 2 is the square of the adiabatic eigen frequency. Here: 


j2 _ (3<r 1 >-4)/ r fl Pd(r 3 ) 

± f* r 2 pd(r 3 ) 

with 


(i) 


<r x 


Pd(r 3 ) 

fr R P d(r3) 


( 2 ) 


The lower bound r of the integration was placed at the bottom of the outer 
envelope, which is in all his cases deep inside the star, mostly very close in 
distance to the centre. Formally, it should be placed at the very centre, at r 
= 0, but truncation to a small r-value is allowed because at the base of the 
envelope the relative amplitude A r/r is already many powers of ten smaller 
than its value at the surface and therefore the deepest regions do not contribute 
to the stellar (in) stability. 

Stothers (1999) gave various examples of the behaviour of stellar models for 
different values of < Tj >. We quote one of them: his Fig. 2 shows the time- 
dependent distance from the stellar centre of various layers of a dynamically 
unstable and pulsationally stable supergiant (log L/L© = 6; T e jf = 10 000). In 
that case a value < T\ > = 1.330 was found, which implies, with eq. (1), a 
negative value of cr 2 . Hence, the model should be dynamically unstable, and 
this is confirmed by the model calculations. An additional result is that the 
model appears to be pulsationally stable. 

Two comments are here in order. First, it appears, from a look at the model 
calculations given, that stars where < T 1 > < or > 4/3 are often also (but in 
rough approximation only) pulsationally unstable or stable, respectively. Next, 
one may wonder whether averaging T\ over the atmosphere only c an also give 
information on the instability of the whole star, because for the stars situated 
in the upper part of the HR diagram the atmosphere occupies a considerable 
fraction of the star. We refer to our Fig. 4 which shows that A r/R (where Ar 
is the depth range between r R = 0.001 and 10) attains values of the order of 
0.1 to more than 0.3. A comparison with Stothers's model calculations shows 
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Figure 1: Areas in the upper part of the HR diagram where stars are dynami- 
cally unstable according to the criterion < T 1 > < 4/3.The diagram shows the 
’yellow-red’ and the ’blue 1 dynamic instability regions. No model calculations 
are so far available above and below the upper and lower dashed horizontal lines, 
respectively. 


that this extent is not large enough, because A r/R is still considerable at tr = 
10 . 

The question that is relevant to the present study is that of the delineation 
of the areas in the HR diagram where < IT > < 4/3. The data for the yellow- 
red hypergiant region are given in Fig. 2 of Stothers and Chin (2001). In 
that paper, particular attention was given to the high-temperature border line. 
Its position, depending on the assumed value for the parameter for convective 
mixing, is given in our Fig. 1. The positions of the lower and upper boundary 
of the yellow-red instability area are tentative. 

The definition of the blue area is more uncertain. It can be read from 
Stothers and Chin (1996), where the instability area, for a metal ratio Z = 0.03, 
is displayed in their Figs. 1 and 2. The red edge parameters are listed in Table 
1 of the quoted paper. The blue area has a high-1 uminoeity cut-off around 
log (L/Lq) = 6.1. This value is uncertain; it may be as low as 5.9. The low- 
luminosity cut-off is around k>g(L/L©) = 5.4. The tentative border lines of the 
blue area are also drawn in our Figure 1. 

Two aspects are brought forward by Figure 1: the yellow-red area contains 
all red and yellow hyper giants. The blue area contains virtually all S Dor stars 
(= LBV’s) in their stable states. This is one of the bases for the assumption 
(Stothers and Chin, 2001, in line with Lamers et al., 1983) that yellow hyper- 
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giants are dynamically unstable stars that are evolving blueward and that, hav- 
ing entered the blue instability region, show up as S Dor stars (LBV’s). This 
assumption implies that S Dor stars are very evolved, dynamically unstable, 
stars. 

3 Photospheric instability defined by g e f / < 0; 
the case of HR8752 

In order to study photospheric instability we have calculated photospheric mod- 
els and derived from these models the four components of the acceleration: the 
Newtonian acceleration, g^\ the radiative value, g r \ the turbulent acceleration, 
Qt\ and the acceleration due to the wind, g w . They add up to g^ff- For stars 
with T e fj below to roughly around 10 000 K, the turbulent and Newtonian 
components are the most important ones; for hotter stars g r contributes too, 
and for the hottest and most luminous ones g w becomes a major component. 
Cases can occur in which their sum, <7 C //, is negative, hence directed outward. 
For cool stars, g n and g t are the competing components. Actually, for the most 
luminous cool objects the absolute value of g t can exceed that of gs , and in 
that case g e ff will be directed outward. We refer to Table 1, further below. 

For the computations, the input parameters of the models are, apart from 
L/L 0 and Tefj- 

- The stellar mass, M/Mq. This value is derived from Maeder and Meynet 
(1988) We took their masses for blueward evolving stars. An interpolation 
programme was developed to obtain M/M 0 for input L/Lq and T c /f values. 

- The stellar rate of mass loss, M. These values were taken from de Jager, 
Nieuwenhuijzen and van der Hucht (1988). We realise that some authors advo- 
cate the use of higher values of M for the most luminous stars (factors 2 and 
even 3 are sometimes mentioned) but we decided to stick to the published data 
for the sake of consistency and also because there are no compelling observa- 
tional reasons yet that suggest higher M values. Anyway, the value of g w , which 
depends on M, is in no case a decisive contributor to g e ff . 

- The microturbulent velocity component in the line-of-sight, Cm* This quan- 
tity is unknown over the greater part of the HR diagram. Since input data 
are needed for a consistent photospheric model, we had a look at the litera- 
ture and compared the observationally derived Cm values with the photospheric 
sound velocity v s at tr — 0.67. We found for the ratio Q^/v a values ranging 
between 1.4 and 2.2, and clustering around 2. Therefore we decided to make 
the computations for two cases: Cm = v a (in order to have a ’minimum case’) 
and Cm = 2t which may be more than a maximum case for the most lumi- 
nous hot stars, because for such stars, for which the sonic point is situated in 
photospheric regions (Fig. 3), the large observed Rvalue is partly due to the 
strong u(r) gradient over the region of line formation, and the real value of Cm 
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is therefore smaller. An example is the B2 supergiant HR80077 (Carpay et al., 
1989) for which = 23 km s~ 1 . For cooler and less luminous stars ( ^ = 2v 8 is 
a valid approximation. 

The various accelerations were derived on the basis of the equation of con- 
servation of momentum, written as 

— dP, GM dM 2 (C„) 2 dlnp 

pdz r 2 r 2dz 2dz 

in which we write 

Id P T 

9r - 

p d z 

where 

Ft = \j Iu) COS 2 

Further 

_ -M 
Vw 4nr 2 p 

Eq. (3) is written as 

9eff -9N +9r+9w + 9t' ( 7 ) 

which defines the four components and their sum. All five quantities are 
depth dependent, and not to a small degree: Cases are rare in which g e /j varies 
by less than a factor 2 over the photosphere and there are many cases in which 
the range is larger than a factor 10. Evidently, the structure of a photosphere 
in which g c f / varies over such a large range differs greatly from one in which it 
is constant with depth. 

To illustrate the relative importance of the ^-components for various combi- 
nations of Ten and L/L©, we refer to Table 1. 

In view of the strong variability of ffc//(^) the photospheric models were 
derived by an iterative method. With the above input values and with an 
estimated initial value for g e ff (we took (<7c//)m = 9n + 9t at t = 0.67) a 
photospheric model was interpolated in the Kurucz set of models. For that 
model the depth dependent values of g e f / were calculated. 

The next step was a calculation of a model with these depth-dependent 
<? c //(t) values, on the basis of a T(r) relation interpolated in the Kurucz models 
for the given value of T e fj and the average g e f f value. For the new model g K f / (r) 
was derived anew and a new (third approximation) model with depth -dependent 
9rff( T ) was derived; a process that was repeated until convergence was reached, 
in the sense that the new < g e ff{j) > did not differ significantly from that of 
the previous step. 


8du). 


( 4 ) 


( 5 ) 


( 6 ) 
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Table 1: The components of g e f j at an optical depth of Tr = 0. 1 in photospheres 
defined by some combinations of T e f/ and log(L/L©), the latter written in the 
table heading as L. The table also gives the corresponding stellar mass M/Af© 
(written in the heading as M) for blueward evolution, the corresponding rate 
of mass loss —log M (written as M), and the assumed microturbulent velocity 
component . Af, M and were derived as described in the text. The 
microturbulence is in kms -1 and the accelerations in cms -2 . 


T '1I 

L 

M 

M 


9n 

9r 

9t 

9w 

9'!J 

7079 

5.45 

7.91 

4.86 

4.94 

1.68 

-.02 

-.32 

0 

1.34 

7079 

6.05 

34.60 

3.32 

4.81 

1.87 

-.02 

-.37 

0 

1.48 

7943 

5.35 

7.50 

5.17 

15.25 

2.55 

-.13 

-.78 

0.01 

1.65 

7943 

5.7 

13.1 

4.47 

7.48 

2.31 

-.12 

-.54 

0.01 

1.66 

10000 

6.0 

32.7 

4.94 

20.11 

7.55 

-.76 

-4.80 

0.05 

2.04 

15850 

5.7 

19.8 

5.24 

13.80 

35.0 

-3.6 

-.7 

1.3 

32.0 


Table 2: A typical course of the iterations for a model with low T e // . Here, r son 
is the optical depth of the sonic point, < g c ff > is the average value of g c f j over 
the tr interval between 0.007 and 0.75, and < g e ff > ou t is the average over the 
depth region 0.007 to 0.2. The model is for T c // = 7079 K and log(L/L©) = 
5.9; = 4.87 km s' 1 . 


iteration # 

on 

< 9eff > 

^ 9c f f >out 

1 

- 

0.79 

0.79 

2 

.0017 

1.55 

1.46 

3 

.0007 

1.44 

1.29 

4 

.0007 

1.47 

1.34 

5 

.0007 

1.46 

1.33 
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Table 3: As Table 2, but for a higher temperature. T c // = 15850 K and 
log(L/L 0 ) = 5.7; = 13.8 kms* 1 . 


iteration # 

Tson 

< 9eff 

< f >out 

1 

- 

90 

90 

2 

.29 

56 

53 

3 

.004 

22 

11 

4 

.11 

34 

22 

5 

.012 

48 

42 

6 

.007 

38 

29 

7 

.012 

49 

44 

8 

.007 

36 

28 


For low values of T cffy roughly for Ten < 10000 K, convergence was nor- 
mally reached in 3 to 4 iterations, as is shown in the example of Table 2. As 
a rule we stopped after the fourth or fifth iteration. For higher temperatures, 
where all four # -components come into play, successive iterations appear to al- 
ternate around an average value, as is shown in the example of Table 3. We 
checked that this average is close to the ’best 1 value and therefore, starting with 
the fourth iteration, we usually took as input parameters for the 2nth iteration 
the average values of g e ff{ t) between those of the (2n — l)th and the (2n — 2)th 
iterations. This procedure worked well. 

For the resulting model two values for the average of g t j* were derived. The 
first, called < g e ff >, is the average over the depth interval r = 0.007 to 0.75; 
the other, < g e ff > 0 ut was averaged over a more outer region: 0.007 to 0.2. 
There are models for which the latter is negative while the former is not. These 
models thus define a transition region between models for which the average 
g e j f values are positive and negative respectively. That transition region is too 
thin to make it well visible in Fig. 2. 

Note in Fig. 2 the position of the region with < g e ff > < 0. It strikes us 
that the three stars mentioned earlier ( p Cas, HR8752 and IRC+ 10420) are all 
situated close to the line < g c ff > = 0. IRC+10420 even lies in the region of 
negative values. This suggests that the low-temperature border line of the area 
where < g e jf > < 0 is an obstacle for blueward moving supergiants. This ob- 
servation implies that the hypergiant characteristics are (at least partly) related 
to their positions in and near that instability area. Earlier, two of us (de Jager 
and Nieuwenhuijzen, 1997) described the repeated to-and-fro movemements of 
HR8752 along the T e ff axis as Trouncing against the border of the instability 
region’. Fig. 2 specifies the location of the bouncing as the low -T c // border of 
the area where g c ff < 0. The object HR8752 is remarkable in that respect. 
Starting around the years 1983 - 1985, its T e ff value has steadily and dramat- 
ically increased till about the year 1998 (Israelian et., 1999), which implies a 
steady compression, reminiscent of a dynamic instability (Nieuwenhuijzen and 


8 




6 0 


66 

64 

c 62 
3 

« 

^ 60 
J 1 5 . 

56 

54 

5.2 


X 

LBV (quiescence) 

+ 

LBV (eruption) 

• 

galactic hypergiant 

O 

in other galaxies 

• 

supergiant 



» \ * . 


45 44 43 42 


4.0 3.6 

»°fl T off 


37 36 3.5 34 



4 5 4 4 4 3 4 2 4 1 4.0 3 9 3 8 3 7 36 3 5 3 4 


•ogTrt, 


Figure 2: The HR diagram with lines of equal < g e f / > values. The upper 
diagram is for = v 3 ; the other is for twice these values. 
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de Jager, 2000; de Jager et al., 2001). We think that this behaviour can be 
interpreted as a dynamic instability triggered by the decrease of < g e jf > to 
below zero. The observations show that < g e ff > <0 the last time in 1978. 
That period was followed by one of enhanced mass loss (around 1980 - 1982), 
and that event was followed, starting in the period 1983 - 1985, by the long 
period of T c // increase. 

4 The depth of the sonic point 

Along the lines described in the previous section we derived the optical depth 
t B oti where the wind velocity equals that of sound. Fig. 3 gives lines in the 
HR diagram of constant r aon values, again for the two assumed values for C,u 
We draw attention to the position of IRC+ 10420 and the neighbouring hyper- 
giant HD33579 (open circle in the diagrams). The latter star cannot be used 
to check the results of Fig. 3, because it is a redwards evolving star, as follows 
from its chemical abundances (Humphreys et al., 1991), and from its large mass 
(Nieuwenhuijzen and de Jager, 2000). The other object though, IRC+10420, is 
an evolved post-red star (Klochkova et al., 1997, Blocker et al, 1999; Nieuwen- 
huijzen and de Jager, 2000). It is rapidly evolving blueward, though not as 
quickly as HR8752: its spectral type changed from F8 I + in 1973 to mid-A in 
1998. Hence its effective temperature has increased by about 2000 K in these 26 
years. We forward the hypothesis that its instability is due to the combination 
of three causes: The star is located at the boundary of the < g e fj > < 0 area 
and inside the r aon > 0.01 area, and the object is also dynamically unstable (cf. 
Fig 1). These effects may explain the instability and the large surrounding gas 
clouds. Blocker et al. (1999) suggest that a phase of heavy mass loss occurred 
some 60 to 90 years ago. 

One aspect of the large r son values is that the replacement time t rcp i of 
the photosphere is relatively short. We estimate it at H/v s where H is the 
scale height. Inserting the expressions for the two variables and taking fj. and 
7 equal to unity (for order-of- magnitude considerations), one obtains t rep i = 
(?RT e ff) l / 2 /g e ff. For a star with T e ff = 10 000 K and g e ff = 1 we have 
trepi = 10 days. 

A related matter is that of the relative extent of the photospheres. To get 
an impression, we calculated for a number of the model atmospheres the value 
of A r/R, where r is the radial distance between the levels with tr = 0.001 and 
10. The outcome is presented in Fig. 4. It appears that for the three unstable 
yellow hypergiants that were mentioned several times in this paper the relative 
extent of the photosphere is larger than 0.3, a considerable fraction of the star. 
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Figure 3: The HR diagram with lines of equal r son values. The upper diagram 
is for = v a at r R = .67. The lower one is for twice these values. 
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Figure 4; Lines of equal Ar/R-values in the upper HR diagram. Ar is the radial 
distance between the levels with Tr = 0.001 and 10; R is the stellar radius. The 
diagram is for the case = 2v 8 . 

5 Conclusions 

The main results of this study are contained in Figs. 1 to 3. They demonstrate 
the presence of two regions of stellar interior dynamic instability, one in the 
yellow-red and another in the blue part of the HR diagram. In addition, it 
appears that the atmospheres of blueward evolving supergiants become unstable 
when their effective temperature has risen to about 8000 K, for two reasons: The 
average g e ff value becomes negative, and the sonic point is getting situated in 
photospheric regions. 

We have compared these results with the recent life histories of two yellow 
hypergiants, HR8752 and IRC+10420, and we have shown that their behaviour 
can be explained on the basis of the data presented in this paper. 

It is obvious that the two areas in the HR diagram, defined by g e f/ < 0 and 
by r aon > 0.01 differ in shape and position from the ’yellow void’ and the ’blue 
instability region’, as described in section 1. That is because these latter two 
regions were defined according to five criteria. We think that the present study 
does not invaliditate the earlier results, but they do allow more specific studies 
of atmopsheric instability. 
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